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: Abstract 



We investigate the effect of resonant temporal forcing on an anisotropic 



■ system that exhibits a Hopf bifurcation to obhquely travehng waves in the 

absence of this forcing. We find that the forcing can excite various phase- 
locked standing-wave structures: rolls, rectangles and cross rolls. At onset, 
at most one of the two - rolls or rectangles - is stable. The cross rolls can 
arise in a secondary bifurcation and can be stable. Experimentally, they 
would appear as a periodic switching between a structure in which the 'zig'- 
component dominates and one with dominating 'zag'-structure. Since there 
are two symmetry-related states of this kind one may expect disordered struc- 
tures to arise due to the break-up of the pattern into domains. The results 
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are consistent with recent experiments on electro-convection in nematic liquid 
crystals by de la Torre and Rehberg. We also apply the general analysis to 
a model of the behavior near a Lifshitz point, where the angle of obliqueness 
vanishes. This analysis indicates that phase-locked standing rectangles are 
always unstable in this parameter regime. 
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I. INTRODUCTION 



Temporally-periodic forcing can qualitatively change the behavior of systems exhibiting 
spatial or spatio-temporal structures. In systems which form oscillatory structures, strongly 
resonant temporal forcing has a particularly large impact. For example, in one-dimensional 
systems that undergo a Hopf bifurcation to stable traveling waves, forcing can excite stable 
standing waves below the threshold for the traveling waves; these standing waves are phase- 
locked to the forcing |[T|J^,^. This theoretical prediction has been confirmed experimentally 
in Taylor-Dean flow and in electro-convection of nematic liquid crystals (EHC) [0. 

In the nematic phase, the liquid crystal system exhibits an axial anisotropy due to a 
partial alignment of the molecules; this defines the director. In the first EHC experiments the 
waves always propagated along the director [^]. Recently the experimental parameter regime 
has been extended to include a regime in which disordered patches of waves, traveling at an 
oblique angle to the director, are observed 0. In this case, temporally-periodic forcing excites 
standing waves in the form of either oblique rolls or rectangles. In addition, more complicated 
structures are found 0. These observations provide the experimental motivation for our 
theoretical investigation. 

This paper examines the effect of resonant forcing on oblique traveling waves produced 
in a Hopf bifurcation of the spatially uniform state of a two-dimensional anisotropic system. 
The symmetry of the system forces coexistence of waves traveling in four directions. The 
possible interactions of these waves makes a complete analysis of the problem extremely 
complicated. However, one of the striking features of the excited roll and rectangular struc- 
tures is their strong spatial coherence. This suggests that as a first step we focus on the case 
where the amplitudes of the different modes are space- independent. We therefore investigate 
an extension of the four complex amplitude equations analyzed previously 0; this model 
takes the temporal forcing into account 0, but neglects large scale spatial variations in the 
amplitudes. This analysis elucidates the connection between the different structures excited 
in the experiments and allows us to make some predictions for future experiments. We also 
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consider the transition from oblique to normal traveling waves using suitably extended cou- 
pled Ginzburg-Landau equations which are valid close to a Lifshitz point where the angle of 
obliqueness goes to zero. 

Formally, our analysis applies to weakly forced waves arising through a Hopf bifurcation. 
We expect, however, that some aspects may pertain to parametrically excited waves in 
systems that do not undergo a Hopf bifurcation, but which support weakly damped waves 
in the absence of forcing. In particular, this approach has been successfully applied to 
surface waves arising in the Faraday experiment with low- viscosity liquids [|l^ and is related 



to methods used in studies of spin waves in ferromagnets [|ri|]. Of course, when applying 
the present approach to these systems only the phenomena arising in the linearly damped 
regime are of relevance {i.e., for parameter values below the Hopf bifurcation point). 

The organization of the paper is as follows. Section II presents the amplitude equations 
for weakly forced oblique waves. Section HI reviews and interprets results on the transition 
from traveling waves to phase-locked standing waves. The phase-locked structures, as well 
as the transitions between them, are discussed in detail in section IV. Section V contains 
our analysis of the Lifshitz equations. The theoretical results are then compared with the 
experimental ones in section VI. The main results of the paper are summarized in the 
conclusion, section VII. 



II. AMPLITUDE EQUATIONS FOR RESONANTLY FORCED OBLIQUE WAVES 

This paper focuses on Hopf bifurcation in axially anisotropic systems in the case where 
there are four neutrally stable modes at the Hopf bifurcation point. These modes correspond 
to waves that travel at an oblique angle 7 to the axis of anisotropy. Thus, a typical scalar 
field such as the vertical component of the velocity in convection, Vz, is given by 

V,{x, y, t) = (2ie^(''^+™) + z2e'^-'"'+Py'> + z^e-'^'^'^+Py^ + z^e-'^-'^'^+^y^) e'""^' + c.c. + h.o.t., (1) 

where we have suppressed any dependence of Vz on the vertical coordinate z, and where 
ujh is the Hopf frequency. Note that zi and z^ correspond to left- and right-traveling plane 
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waves in the direction of gx + py; similarly Z2 and correspond to oppositely traveling 
plane waves in the direction of — gx + py. Thus, the waves travel at an oblique angle 
7 = tan~^(p/g) to the a;-axis (the axis of anisotropy). Here we assume that the amplitudes 
Zi are complex-valued functions of a slow time T, but that they do not depend on the 
spatial coordinates. The time evolution of the amplitudes is determined by a system of 
coupled first order ordinary differential equations which can, in principle, be derived from 
the governing hydrodynamic equations (with periodic boundary conditions) using center- 
manifold reduction However, in EHC this is not yet possible; even the origin of the 



Hopf bifurcation has not been understood at the hydrodynamic level . Nevertheless, the 
general form of the amplitude equations can be determined from symmetry considerations 
assuming that the waves are due to a Hopf bifurcation of the motionless state as experiments 
suggest [|1^]. Consequently, most of our results are not specific to the EHC-experiments, but 
pertain to a larger class of resonantly forced systems with the same underlying symmetries. 
In a more heuristic fashion, we also expect our results to be relevant to surface waves in 
the Faraday experiment |T^. We will, however, focus on the EHC-system in order to make 
contact with available experiments 0. 

We consider a general system that is invariant under translations in the plane; since our 
solutions are spatially doubly-periodic, the group of translations corresponds to a torus T^. 
Moreover, we assume reflection symmetries in planes parallel and perpendicular to the axis 
of anisotropy; these reflections are denoted by ki and /t2, where ki : {x,y) {x,—y) and 
K2 : {x,y) {—x,y). These spatial symmetries act on the amplitudes as follows: 
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Here the phase shift 61 results from a translation {x,y) — > (x + di/2q,y + 9i/2p) and the 
62 phase shift results from a translation by (—6*2/25, 6'2/2p). The amplitude equations are 
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equivariant with respect to these transformations. In the absence of resonant temporal 
forcing the normal form of the amplitude equations possesses an additional phase shift 
symmetry that we interpret as a time-translation symmetry t —>■ t + (p/ooh- It acts on the 
amplitudes as follows: 

■ z = e'^^z, (3) 

where z = {zi, Z2, z^, Z4). An analysis of the cubic truncation of the amplitude equations, 
equivariant with respect to the symmetries (H) and @, is presented in 

The temporal forcing destroys the phase-shift symmetry (^), thereby allowing additional 
terms in the amplitude equations. The frequency of the external forcing is c^e = + Au, 
where we assume that the detuning Au is small and that n = 1 or 2 (n = 2 in the EHC 
experiments). Moreover, we assume that the temporal forcing is weak and therefore we keep 
only the lowest order term which does not satisfy the time-translation symmetry (j^). The 
amplitude equations truncated at cubic order are 
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= uzi + fiz^ + (a\zi\ 
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There are three small unfolding parameters: the distance from the Hopf bifurcation in the 
absence of forcing, Re{i'} = z/,,, the detuning Auj of the external frequency with respect to 
the Hopf frequency, I'm{u} = z/j oc Auj, and the forcing amplitude F, \fi\^ oc F"^ (In 
the case of parametrically excited, weakly damped waves, is proportional to the linear 
damping of the waves; it is always negative.) Without loss of generality, we choose the 
temporal phase of the forcing such that fi is real and positive; all other coefficients are 
complex. Finally, we replace Uh by Ue/n in the expression ([1|) for the velocity component 
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In the absence of forcing (/i = 0) eqs. (^) possess up to seven branches of periodic 
solutions that bifurcate from the origin at z/^ = 0. Three of these correspond to different 
types of travehng waves and four to (superpositions) of standing waves. In particular, they 
have the form 

/: z = (t;, 0,0,0) TRo, 
II : z = {v,v,0,0) TRe^, 

III : z = {v,0,0,v) TRe\ 

IV ■ z = (t;,0,t;,0) SRo, (5) 
V : z = {v,v,v,v) SRe, 

VI : z = {v,iv,v,iv) ARo, 
VII : z = (f 1, V2, f 1, V2)SCR. 

The abbreviations denote traveling rolls (TRo), two kinds of traveling rectangles distin- 
guished by their direction of propagation relative to the axis of anisotropy {TRe-^,TRe^^), 
standing rolls (SRo), standing rectangles (SRe), alternating rolls (ARo) and standing cross 
rolls (SCR). The standing cross rolls, which are characterized by two complex amplitudes Vi 
and f 2, do not always exist. The numbering is the same as that used in our previous work on 
the unforced system p|. In addition to these periodic solutions one finds more complicated 
(branches of) attractors such as heteroclinic cycles connecting three of the above solutions 
i- 

In the presence of forcing, oppositely traveling waves are linearly coupled (see eq. (^)) 
so that pure traveling waves do not exist. For example, the forcing destroys pure right- 
traveling Ti?e" by coupling it to the left-traveling Ti?e"; thus the Ti?e" are replaced by a 
general superposition of left- and right-traveling rectangles, i.e. by (f i, f2, "^2, f i). Similarly, 
forcing destroys the ARo solution, which is an equal-amplitude superposition of standing 
rolls which differ by tt/2 in their temporal phase. In the next section we see that the pure 
standing waves {SRo and SRe) persist and that they can phase-lock to the temporally- 
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periodic forcing. 



III. TRAVELING VS. PHASE-LOCKED STANDING WAVES 

In this section we study the effects of resonant forcing on competition between vari- 
ous travehng-wave structures and their standing-wave counterparts. Specifically, we focus 
on the stability of the phase-locked standing waves to traveling-wave disturbances. In the 
oblique regime there are two simple types of phase-locked standing waves {PSW)\ phase- 
locked standing rolls {PSRo) and phase- locked standing rectangles {PS Re). The stability 
of the PSRo state to traveling disturbances is determined by restricting to the subspace 
where z = (t>i, 0, f2, 0). Similarly the stability of PSRe to traveling disturbances is deter- 
mined by restricting (H), in turn, to each of the subspaces where z = (fi, f i, f2, "^2) and 
z = (f 1, f2, ^"2, f 1); in the former case the disturbances travel perpendicular to the director, 
whereas in the latter case they travel along it. In each of these three cases, the restricted 
dynamics has the form 

Vl = uvi + /if2 + + /3fi|f2|^, 

V2 = Z/f2 + fivl + l3v2\Vi\'^ + af2|f2|^, (6) 

where 

a = a, (3 = c for Ti?o, (7a) 

a = a + h, j3 = c + d + f for Ti?e^, (7b) 

a = a + d, p = b + c + f iovTReK (7c) 

Here the subspaces are denoted by the traveling waves that are present in each of them when 
there is no resonant forcing (c/. (^)). Eqn. (^ was studied previously in the context of 
periodically forced waves in one dimension [|T],^^; in these papers the phase-locked standing 
waves (PSW), as well as traveling (TW) and unlocked standing waves (SW) were investi- 
gated. The results apply directly to the present system by an appropriate identification of 
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TRo and TRe^^'-^ with the travehng wave TW, and the new phase-locked solutions PSRo 
and PSRe with the standing wave PSW. 

In fig.|I| we sketch a typical phase diagram associated with (H), for fixed detuning Ui, 
in the case that the unforced traveling waves bifurcate supercritically and are stable. In 
the following discussion we employ the notation {TW, PSW) used in the one-dimensional 
case. For weak forcing the trivial equilibrium loses stability via Hopf bifurcation along the 
line marked H; this bifurcation leads to stable traveling waves, TW, for positive values of 
Uj.. (The unlocked standing waves, SW, are unstable in this case.) Along the line marked 
PSW, which is given by /x^ = |z/p, a steady bifurcation of the zero solution of (H) occurs; 
this bifurcation produces the phase- locked standing waves PSW with frequency Ue/n. The 
amplitude r = = \v2\ of the PSW satisfies 

\N\^r^ + 2r^Re{iy*N} + \iy\^ - 1^^ = 0, (8) 

where N = a + p. If Re{u*N} = v^Nr + ViNi > the bifurcation to PSW is forward, 
otherwise backward. Thus, depending on the detuning z/j, the bifurcation to the PSW can 
be in either direction, independent of the bifurcation direction of the standing waves in the 
unforced Hopf bifurcation. The saddle-node line of the PSW, marked SN in fig.|l], is given 

by 

fi' = {Im{u*N}y/\N\^ (9) 

for Re{u*N} < 0. 

The linear stability analysis of the PSW is simplified by the fact that the perturbations 
(^1,^2) of the steady solution (^1,^2) = {v,v) of (H) fall into two classes: Ui = U2 and 
ui = —U2. We focus first on the perturbations that preserve the standing- wave character of 
the waves {ui = U2). The determinant and trace of the resulting 2x2 stability matrix are 

Det = \u + 2Nr^\^ - (10) 
Tr = 2{ur + 2Ny), (11) 
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where solves (H). The determinant vanishes at the saddle- node bifurcation (^). There is 
a secondary Hopf bifurcation when the trace vanishes, provided the determinant is positive. 
This bifurcation produces unlocked standing waves SW. The perturbations of the form ui = 
—U2 include the marginal translation mode and a parity-breaking instability to traveling 
waves. The parity-breaking instability occurs when the trace of the associated Jacobian 
matrix vanishes, i.e. when 

Tr = 2{ur + 2ay) = 0. (12) 

It is denoted by PB in fig.0. A typical bifurcation diagram is presented in fig.^ As long 
as the unforced traveling waves are stable with respect to the standing waves {Pr < 0) the 
secondary Hopf bifurcation of the PSW to unlocked standing waves is always preempted by 
the parity-breaking bifurcation to traveling waves. In addition, modulated waves exhibiting 
three different frequencies are possible; this follows from an analysis of eq. (P) in the vicinity 
of the point where the Hopf bifurcation and the steady bifurcation merge (/x^ = uf, = 0) 
||T|]. This CO dimension-two point corresponds to a Takens-Bogdanov bifurcation with 0(2)- 
symmetry, which was analyzed by Dangelmayr and Knobloch 

Thus the stability of the phase-locked states {PSRo, PSRe) to traveling disturbances 
{TRo, T-Re-*-' II) and to unlocked standing waves (SRo, SRe) is understood in some detail 
based on previous studies. The competition between the different kinds of phase-locked 
standing waves is the subject of the next section. 
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IV. COMPETITION BETWEEN PHASE-LOCKED STANDING ROLLS AND 

STANDING RECTANGLES 



The relative stability of PSRo and PSRe is particularly relevant to our understand- 
ing of the EHC experiments since both states are observed. This stability question can 
be addressed by restricting (Q) to the standing-cross-rolls subspace where {zi, Z2, z^, Z4) = 
(f 1, f2, f 1, t>2)- Note that this subspace contains both the PSRo {vi = or = 0) and the 
PSRe {vi = V2) solutions. The amplitudes vi and V2 satisfy 

vi = vvi + fivl + {a + c)vi\vi\'^ + {b + d)vi\v2\'^ + fvlvl, 

V2 = VV2 + n,vl + (6 + d)v2\vi\^ + {a + c)v2\v2\^ + fvlvl- (13) 

These equations differ from those describing the forced 0(2)-Hopf bifurcation (H); the modes 
with amplitudes vi and V2 are uncoupled at the linear level since they represent standing 
waves with different orientations. Moreover, the nonlinear coupling term with coefficient / 
is absent in (H). 

We focus first on the stability of the PSRo solution, which has the form (f 1, ^2) = (m, 0), 
to perturbations in the (0, f2)-direction. The resulting 2x2 stability matrix DpsRo has the 
following trace and determinant: 

Tr{DpsRo) = 2{Ur + {br + dr) (14) 

Det{DpsRo) = W+{b + dylsRol^ - W + {a + c - fylsRol^- (15) 



The amplitude trsro satisfies (H) with N = a + c. Provided the determinant (|15D is positive, 
there is a Hopf bifurcation at = —{br + dr)rpsRo't thus the presence of the PSRo in one 
direction leads to a shift of the primary Hopf bifurcation in the other direction. A steady 
state bifurcation occurs for 

rlsRo _ Re {iy*{a + c - b - d - f)} 

2 ^ \b + d\^ -\a + c- ■ ^ ^ 

We will show that this bifurcation leads to phase-locked standing cross rolls PSCR which 
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are also pertinent to the EHC experiments. The PSRo are stable in the (0, f2)-direction as 
long as the trace (0) is negative and 

2Re{v*{a + c-b-d- f)} < {\b + dp - |a + c - /p) rlsRo- (17) 

Note that p depends on both the control parameters and Vi. Moreover, since the amplitude 
rpsRo can be varied independent of u by changing the forcing amplitude /i, eq. ([l7D reduces 
at onset to 

2Re{iy*{a + c-b-d- f)} <0. (18) 

Finally we note that the symmetry of the PSRo solution results in a doubling of the eigen- 
values associated with perturbations at an oblique angle to the rolls (see, for example, 0). 
Hence the stability calculations of this section actually determine four of the eight eigenval- 
ues of the full stability matrix associated with (^) linearized about PSRo; the remaining 
four eigenvalues were determined in the previous section. 

We now examine the stability of the phase-locked standing rectangles (PSRe) in the 
standing-cross-rolls subspace. The PSRe solution is an equilibrium solution of (|T3p of the 
form (fi,f2) = {u,u). We examine its stability with respect to perturbations of the form 
{w, —w) by substituting (t>i, V2) = {u + w,u — w) in ([13|). The resulting 2x2 stability matrix 
DpsRe has the following trace and determinant: 

Tr{DpsRe) = 2Ur + 4{ar + Cr- frYlsRe. (19) 

Det{DpsRe) = |z/ + 2(a + c - fYlsRe? - ^ + 2(6 + dyl^ReW (20) 
where rpsRe is a solution of (|]) with N = a + h + c + d + f. Note that 

Det {DpsRoirlsRo)) = -Det {DpsRe{2rlsRe)) ■ (21) 
Consequently a steady bifurcation to PSCR occurs at 

rlsRe = P, (22) 
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where p is defined in ([T6|). The coincidence of the steady state bifurcations of the PSRe— 
and the P5'i?o— branches, both to PSCR solutions, suggests that they may be connected by 
the same branch of PSCR; this provides a mechanism for a continuous transition between 
PSRo and PSRe. We will show, in a particular limit, that this possibility can occur. 

The PSCR are steady state solutions of (|TE|) with V1V2 7^ and It'll 7^ |f2|. We let 
f 1 = i?e** and f 2 = re**^, and eliminate the phases from (|l^) to obtain (for r 7^ 0) 

„+(^a + cy + {b + d)R' + fR^^-^—l±^ 23 

z/ + (0 + d)r^ + (a + c — j)R^ 

Moreover, if r and R are both non-zero, then 

\u+{b + d)R^ + (a + c - /)r2|2 = |iy + (a + c - + (6 + c/)r2|2. (24) 

This has a solution = R"^, which corresponds to the PSRe, and also a solution 

r^ + R^ = 2p. (25) 

Substituting = 2p — i?^ in (^) leads to a quartic equation for i?^: 

CsR^ + CqR^ + CiR"^ + C2R^ + Co = 0. (26) 

The coefficients Cq, C2 and Cg are given by 

Co = |z/ + 2(6 + rf)p|2 (|z/ + 2(a + c)p|2 - /i^) , (27) 
C2 = 2Re [{p + 2(6 + d)p) {h* - /*) (2p(z/ + 2(a + c)p)(/i* + /*) - p') } , (28) 
C8 = |/^'-/T, (29) 

where h = a + c — b — d. We use the reflection symmetry that interchanges R^ = |fip and 
= 1^2 P to simplify this equation. In particular, we observe that for each solution Rf there 
must also exist a solution R2 = 2p — R\ {cf. (|25|)). We let 

u = R^{2p- R^) (30) 

and write the solutions to (^) as 
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-Ch±JCi-16CaC,p' 



which yields the symmetry-related pair of solutions 



R^ = p±^p'^- u. (32) 

Note that both u and E? must be positive for PSCR to exist. From (PH), it follows that 
there are at most two distinct PSCi^-solutions of (0); these annihilate in a saddle-node 
bifurcation when C| = IQC^C^p^. We note that the PSCi^-solution merges with the PSRo 
when M = and with PS Re when u = p^. We now consider two limiting cases in which we 
can show explicitly that the PSCR secondary solution branch that connects the 

PSRo- and PS* i?e-solutions. 

The phase-locked standing waves arise through a steady bifurcation at |z/p = /x^, when- 
ever Ur < 0, i.e., for values of Ur below the unforced Hopf bifurcation point. We can use 



center manifold reduction ||T2[ to reduce the complex equations ([T3|) to two real equations 
for the standing- wave amplitudes [1^. For |z/p = /i^, z/^ < 0, the neutral eigenspace of 



(p!3|) linearized about the origin is spanned by Vj = rji = + u*), j = 1,2; the damped 
eigendirections associated with the double eigenvalue 2h'r are given by Vj = 772 = i^ji — u). We 
introduce the slow time r = e^T and let Vj = erjiAj + e^{riiAj3 + ri2Bj3) + h.o.t, p, = |z/| + e^/i2- 



At cubic order in e we obtain 

aAi + aAl + pAiAj, 



dr 

^ = aA2 + aAl + PA2AI (33) 
dr 

where 

II I 1 2 I 1 2 

~^\!nLRe{u*{a + c)}, P=^-^Re{iy*{b + d + f)}. (34) 

Equations of the form (^) also apply to systems with large damping whenever sufficiently 
large forcing is applied. In particular, these equations describe the oblique standing-wave 
components that have the correct temporal phase, relative to the periodic driver, to be 
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excited by the forcing. In contrast, eqs.(^ become invalid for large damping. Note that in 
general the amplitudes in (^) are complex due to translation symmetry; the specific choice 
for the standing-cross-rolls subspace {zi = -23,-22 = 2:4) that leads to eg. (p!3|) breaks that 
symmetry. 

Equations of the same form as (^) arise in the analysis of the 0(2)-Hopf bifurcation 
problem and have been studied extensively [|l^. Generically, the only steady solutions of 
(^) have either equal amplitudes {Ai = A2), or one of the amplitudes vanishes. Thus, at 
onset one has only PSRo (one vanishing amplitude) and PS Re {Ai = ^2). Only in special 
degenerate cases do we expect more complicated small amplitude equilibrium solutions near 
onset. In particular, if a ^ (3, then higher order terms must be included in the normal 
form (|33D and there can exist a small-amplitude mixed-mode equilibrium solution with two 
different non- vanishing amplitudes ]19|,^. The mixed-mode solution corresponds to the 
PSCR. Note that p oc {a — (3) so that p is small in this degenerate situation. This is 
consistent with (^) and ( p2D which show that the secondary bifurcations to PSCR occur 
when the amplitudes t'^psro ^^^d rp^j^^ are of order p, respectively. In the two-dimensional 
control parameter space (/i2, a — (3) the regime in which the PSCR exist is delimited by the 
lines at which the PSRo and PS Re become unstable to PSCR. The stability of the PSCR 
solution depends on higher order terms which were neglected in deriving the amplitude 
equations (D 0,10]. 

The center manifold equations (|33|) lose their validity when i/^ — since the eigendi- 
rections with eigenvalue 2i>r can no longer be eliminated adiabatically; these standing wave 
modes become independent dynamical variables. For nonzero detuning i/j this regime can 
be analyzed by considering the Takens-Bogdanov double-zero point (z/^ = 0,/i = po = |z/j|) 
at which the Hopf frequency goes to zero. Near this point simplification is not achieved by 
a reduction in the number of degrees of freedom. Instead, the form of the linear operator 
allows certain nonlinear terms to be removed p]2| , |2T| . We begin by expanding the amplitudes 



along the two directions Ci = (■s + i)PQ and C2 = 1 in the complex plane, where s = sgniyi) 
distinguishes between the cases of positive and negative detuning. Let 
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Vj = n{CiAj + C24)' J = 1. 2, < n < 1. (35) 

The new dynamical variables Aj and I3j are real. At the double-zero point the linear part 
of flTSD is 



Aj = sBj, (36a) 
4=0. (36b) 

Since s 7^ 0, we can perform a near-identity nonlinear coordinate transformation that re- 
moves most of the cubic terms in the normal form. In particular, we let 

Ai = Ai+ r?{cxA\ + c^BxM + cgis + c^B^M + c^DK 

Wic^Bl + cr^iii + c^A\ + c^B2A2 + cx^Bl)Bx, (37a) 

fii = fii + r?[diA\ + c/a^iii + rfgi^ + d^B^h + 4^^)^! 

W{deBl + c/y^iii + + f^g^^ + dioB^)Bi, (37b) 



with and B2 defined by interchanging the 1 and 2 subscripts in (|37a|) and (|37b| ), re- 
spectively. We choose the coefficients cj and dj so as to remove the cubic terms in the 
Aj equations and four of the ten cubic terms in the Bj equations. Moreover, in order to 
express the unfolding parameters of the resulting equations in terms of the original control 
parameters, we expand z/ and fi as 

/i = /io + e/ii, u = isfio + '^{Ki + ii^ii), (38) 

where < e <C 1. The transformed equations are 

Ai = sBi + e(z/,.i + svii - fii)Ai + eunBi/fio + C(en^,n^), (39a) 
Bi = e^Uri - sua + fii)Bi - 2e/io(z/ii - sfii)Ai 

+4n2{-/i|^(ai + Ci)Al + 2/2^(0, + Cr)BiAl - /i|^(6i + di + fi)AiAl 

+fll{br + dr + 2fr)B2AiA2 + nl{K + dr)BxAl + /io(6i + di- 2fi)B2BiA2} 

+0{en\n'^), (39b) 
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with similar equations for A2 and -82- We focus on the stationary solutions of (|39a|J39bD 



which correspond to the phase-locked standing waves observed in the experiments. In this 
case, ( p9a| ) yields Bj = 0{e,n^), j = 1,2, which takes all terms involving Bj in ( |39b[ ) to 
higher order. We obtain the following algebraic equations for the Aj-. 

= -2/io[e(z/,i - sfii) + 272^0 + Ci)Al + {hi + di + fi)Afj]Ai + 0(e^ en^ n^), 

= -2/io[e(z/a - s/ii) + 2nVo + Ci)^2 + + d, + U)A'^]A2 + 0(e^ en^ n^). (40) 



Note that the equations for the equilibrium solutions of ( p^ ) are of the same form as (^01). 
Thus, as before, PSRo and PSRe are the only phase-locked standing waves in the vicinity 
of the Takens-Bogdanov point, unless -|- q ~ bi + di + fi, in which case higher order 
terms must be retained. In this case, it is known that the existence region of PSCR lies 



between steady state instability lines of PSRo and PSRe pO[- Again this degenerate case 
corresponds to p ~ 0. (Recall that z/,. ^ near the Takens-Bogdanov point.) 

We can combine the above results to obtain the following picture. The steady state 
instability at /i = |z/| produces both PSRo and PSRe, with at most one of the two states 
being stable. This follows from (^TJ); a positive determinant is necessary for stability and for 
''"psRoPSRe ~^ the determinants of DpsRo and DpsRe have opposite signs. (Recall that 
can be decreased by varying the control parameter p, while holding v fixed.) In addition, for 
sufficiently negative, i.e. below the onset of unforced traveling waves, the trace of DpsRo 
(|1^) and DpsRe ([T9| ) are both negative. For p sufficiently small, the steady bifurcations 
from PSRo and PSRe to PSCR can occur before the Hopf bifurcations associated with 
Tr{DpsRo) = and Tr{DpsRe) = 0. The difference in the steady state bifurcation values 
for each of the solutions is given by 

A 2 _ 2 2 
Zi/i _ PpsRo^PSCR f^PSRe^PSCR 

= p2(2(|6 + d|2- |a + c-/|2) - (\a + c + b + d + f\^ -A\a + c\^)) . (41) 



The two bifurcations coincide for p = in which case they lie on the neutral curve p^ = 
For p < there is no bifurcation to PSCR. For p > there are two main cases distinguished 
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by the sign of|6 + (ip — |a + c — / p (c/. (|T^). In these two cases, the conditions for stabihty 
of PSRo and PS Re, respectively, to PSCR are given by 

I + |a + c-/p > 0) : PSRo: rlsRo>'^P, PSRe : rlgj^^ < p, (42a) 

II + |a + c-/|2 < 0) : PSRo: rlsRo<'^P, PSRe : rlgj^^ > p. (42b) 

Figures]^ and § present possible bifurcation diagrams when p > 0, and when both PSRo 
and PSRe bifurcate supercritically. In case I, the PSRe are stable at onset, but become 
unstable to PSCR at larger amplitude {i.e. larger values of /i^). For Ap^ > 0, the PSCR 
merge with the PSRo for even larger /i^ and stabilize them in turn. Our analysis does not 
determine whether the pitchfork bifurcation to the PSCR is supercritical or subcritical. 
In the simplest scenario, the bifurcation is supercritical when A/z^ > and the PSCR 
solution is stable in the SCR subspace. For Ap"^ < 0, the transition from PSRe to PSRo 
is hysteretic and (in the simplest case) the PSCR are unstable. In case II the PSRo and 
PSRe interchange their roles and stable PSCR are expected for A/i^ < 0. If the PSRo or 
PSRe bifurcate subcritically from the basic state more complicated bifurcation diagrams 
are possible. 

We emphasize that p depends on the control parameter u. Within the present framework, 
one therefore expects that it may be possible to choose the detuning z/j such that p is in a 
suitable range to observe the transition to PSCR below the onset of unforced waves (z/^ < 0) 
and before other instabilities {e.g., parity breaking) set in. In particular, by choosing z/j such 
that p is small and positive it may be possible to observe the transition between PSRo and 
PSRe by increasing p. Whether this transition is made via a stable or unstable branch 
of PSCR solutions will depend on the nonlinear coefficients (see (0)). Similarly, if both 
PSRo and PSRe bifurcate supercritically and Ur is sufficiently negative, then p determines 
whether PSRo or PSRe are stable at onset. In particular, PSRo arise stably at onset if 
p < in case I and if p > in case II (see (0) and {^^). As will be discussed below, the 
range over which the detuning and therefore p can be varied effectively is, however, limited 
in spatially extended systems. 
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V. TEMPORAL FORCING IN THE VICINITY OF A LIFSHITZ POINT 



Experiments on the liquid-crystal system have identified a parameter regime where rolls 
travel along the axis of anisotropy as well as a parameter regime where (disordered) waves 
travel at an oblique angle to it. This is reminiscent of the situation of steady normal and 
oblique rolls for which the transition between the two regimes is usually continuous and the 
angle of obliqueness, 7 = tan^^(p/g) (c/. (|^)), goes to zero at the Lifshitz point in parameter 
space p3|j2^ . In the present case of Hopf bifurcation to traveling rolls, the disordered 
character of the unforced waves has made experimental confirmation of the existence of such 
a Lifshitz point problematic; no such experimental verification exists to date. Nevertheless, 
such a point is of special interest since it provides a natural mechanism for a transition 
between the normal and oblique regimes. At the Lifshitz point, the critical modes with 
amplitudes Zi and 2:4 in (|l|) become indistinguishable as ±p 0, as do the critical modes 
with amplitudes Z2, z^; the description in terms of four ordinary differential equations is no 
longer appropriate. Instead, the behavior in the vicinity of the Lifshitz point is modeled by 
a generalized Ginzburg-Landau equation in which the small wave number p is incorporated 
into a slow variable in the y-direction. This description of the Hopf bifurcation in the vicinity 
of the Lifshitz point is analogous to the model of the Lifshitz point in the steady state case 



25| . The latter leads to a single partial differential equation (in two space dimensions) with 
real coefficients, whereas in the Hopf case the system is described by two coupled complex 
partial differential equations. 

For a Hopf bifurcation with translation and reflection symmetry two complex amplitudes 
are required to model the behavior of both the left- and right-traveling waves. We replace 
(I) by 

V,{x, y, t) = e^A{x, y, t)e^K*7w) + ^^b^x, y, t)e^("«*/-P) + h.o.t. + c.c, (43) 

< e ^ 1, where x = e'^x, y = ey and t = eH are slow variables. The spatial symmetries 
of the system result in the following equivariance properties of the generalized complex 
Ginzburg-Landau equations: 
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x-^x + e/q: (A, S, 9,, dy) ^ {Ae'\ Be-'', d„ dy), (44a) 
x^-x: {A,B,d,,dy) ^ iB,A,-d,,dy), (44b) 
y ^ -y : {A, B, d,, dy) ^ (A 5, 9,, -9,). (44c) 

One then obtains the following equations (through order e^) 

dtA - vd^A = ddlA + sd^A + gd^dlA + hd^A + vA + jjB* 

+cA{\A\^ + \B\^)+gA\B\^, 

dtB + vd^B = ddlB + sd^B - gd^dlB + + + ^lA* 

+cB{\A\'' + \B\'') + ~gB\M\ (45) 

where, without loss of generality, we assume that the group velocity v and the (scaled) forcing 
amplitude /i are real; all other coefficients are complex. The Lifshitz equations ( ^5] ) without 
the temporal forcing {i.e. /i = 0) were introduced in [^j and partially analyzed in |p. In the 
unforced case the Lifshitz point is given by Sr = Vr = 0, with the unsealed deviation from 
the Lifshitz point given by (e^s^, e^z/,,). As before we assume the temporal driving is small 
and therefore keep only the lowest order terms /i-B* and ^A* that destroy the phase-shift 
symmetry {A,B) — e'^^{A,B). Specifically, the unsealed forcing is /I = e^/i. Similarly, the 
complex coefficient v has been scaled by e^, i.e. v = e^v. It should be emphasized that all 
terms in (^) are not generally 0{\) in e. In particular, if the (unsealed) group velocity (w) 
and the unsealed imaginary part of the s (sj) are 0{\\ then v and Sj are (9(l/e^). Formally, 
we assume that v and Sj are Oie^) so that v = e^f , s = e^s. In an asymptotic analysis of the 
general case, where v and Sj are 0{1), one would introduce a faster time ti = eH and first 
solve the O(e^) equations before introducing those at 0{e^) p^j27|] . We expect to capture 
part of this general case by assuming v and si are large {i.e. (9(l/e^)) in the final results. 

Without temporal forcing, the growth rate a of modes A,B(x e^**^'''^"'"™^ is given by 
a = Vr — drq^ — SrP^ + giqp^ + /i^p^. Thus, for > 0, a Hopf bifurcation to normal traveling 
rolls occurs at 

Vr = 0, qc = 0, Pc = 0. (46) 
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For Sr < the bifurcation to oblique traveling waves is given by 

1 2 9i^r 2 2drSr „ / a'7\ 



Note that stability of the trivial state with respect to large wave numbers q and p requires 
that dr > 0,hr < 0. Moreover, we assume Ahrdr + gf < so that the instability to oblique 
waves occurs when Sr < 0, with z/,. < in (^). 

With resonant temporal forcing, the location of the Hopf bifurcation to normal, as well 
as to obhque, waves remains unchanged for small /i, i.e. for < < \i^eff{qc,Pc)\, where 

Veffiq^p) = + iqv — dq^ — sp^ — igqp^ + hp'^. (48) 



The frequency on the Hopf bifurcation surface is given by = y |z^e//p ~ /^^- Recall that this 
frequency is the difference between half the driving frequency (for n = 2) and the frequency 
of the critical mode. Note that the frequency uo decreases as /i increases, until it goes to 
zero for /x^ = ii^{q,p) = |i^e//P- At such points, the Hopf bifurcation is replaced by a steady 
bifurcation of the trivial solution of the amplitude equations (^Sj). This steady bifurcation 
leads to (standing) waves which are phase-locked to the forcing. 

To get some insight into the nonlinear solutions of ( ^Sf ) we consider two limiting cases. 
Since the phase-locked waves arise through a bifurcation involving a real eigenvalue, the 
coupled complex equations (|45|) can be reduced to a single equation sufficiently close to the 



neutral surface IT^]. Depending on the kind of extremum (qcPc) chosen as expansion point, 
different equations are obtained. Here we concentrate on a situation that eventually leads to 
an equation of the same form as for the steady bifurcation at a Lifshitz point pSf . Another 
limiting case is obtained by considering oblique waves within (^), i.e. assuming \p\ > 
for small u and fx. This allows us to derive the amplitude equations (^ from the Lifshitz 
equations. 

We first describe the reduction to a single real Lifshitz equation. The condition for a 
Lifshitz point for the standing waves on the normal-roll branch is 

5pV(gc,Pc = o) = o (49) 
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where Qc is determined by dqfj,'^{qc,Pc = 0) = 0. These two conditions determine the critical 
wavenumber q = Qc and the Lifshitz point in parameter space Sr = sf'^, respectively: 



-2{ur - drql)drqc + {vqc + z/j - diql){v - 2diqc) = 0, 
LP {ui + vqc - diql){si + Qrqc) + giqMrql - Vr) 



(50) 
(51) 



^ drql — Vr 

Note that the presence of the forcing has shifted the Lifshitz point for the standing waves 
to a non-zero value of s^- This is in contrast to the unforced case where the Lifshitz point 
occurs at = with qc = Pc = 0. For large values of v and Si, of order e"^, this leads to 

2h>rdr — l^idj 



+ 0(e«) 



(52) 



and a shift of the Lifshitz point to 



s^/ = --Ug,-2dr'^]+0{e% 

V \ V ' 



(53) 



Thus, in this limit of large group velocity v in p5|), the steady state neutral surface fi^ = fi^ 
has only a single minimum with p = 0; this occurs at fi"^ = fi^ = fiQ{qc,Pc = 0), where qc is 
given by (|52D. (In the special case that v is not large, this need not be the case as will be 
discussed below.) 

To derive the reduced Lifshitz equation from eg. (|i5|) , we introduce super-slow space and 
time scales and expand /i about its critical value /ic = |z/e//(Q'c,P = 0,s^^)|, Sr about its 
critical value s^^, while holding z/ fixed (z/,. < 0). In particular, we let 



fi = fic + 5^/i2, Sr = s^^ + 6sri, X = 6x, Y = 6^^"^];, T = 6H, 



(54) 



where < 6 <^ 1. We use the two eigenvectors (/i, — i/g//)* and (/i, z/*jj)* of the linear 
operator evaluated on the neutral surface fi^ = [z^e/zP, and expand the amplitudes A and B 
as follows: 



^a] / , ^ 

= 6 



f ' \ 



\ ^eff ) 



526^^=^ ^ +0(5^). (55) 



At 0[p) we obtain 
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B2 = {{v + 2iqJ)dx + (s^^ + iqcg)dl} A^, (56) 

where i'eff,r = R^i^eff} and Veff is evaluated at (g,p) = (q'cO), and = s^-'^ + isi. Pro- 
jecting the C(5^)-equations onto the left-zero eigenvector (/i, —Ueff) the Lifshitz equation, 

ar^i = {Kidl + ^2^2 + iK^dxdl + K49^)Ai + SAi + TA^\Ai\\ (57) 

is obtained, where the coefficients {j = 1, ...,4), S and F are all real. We find that 



2 



r = -^i?e{.:,,(2g + ^)}. (58) 
The coefficients Ki are given by the appropriate derivatives, with respect to the wave numbers 



q and p, of the growth rate a = z/e//,r + -y/yU^ — i^eff,i^ ^-9- = d^a^qcPc = 0). Note that the 
direction of bifurcation of the spatially uniform, (nontrivial) steady solution of (|57D depends 
on Veffi as well as the original nonlinear coefficients c and g. 

As mentioned before, eq.(|57|) has been studied in the case where the bifurcation to the 



steady uniform solution is supercritical {i.e. in the case where F/S < 0) [^,^. These 
investigations determined that normal and oblique rolls can be stable in two-dimensional 
wave-number regions, where the shape of these regions depends on the parameter values. 
In the present context, such solutions correspond to normal phase-locked standing rolls 
(NPSRo) and oblique phase-locked standing rolls (OPSRo). In the oblique-roll regime 
stable zig-zags arise which connect OPSRo of opposite angles with respect to the director. 
In addition, undulated solutions can also be stable in small islands in wave-number space 



25| , |28|| . In the subcritical case, suitable nonlinear gradient terms in X- as well as ^-direction 
would have to be included when going to fifth order. They would allow the direction of 
bifurcation of NPSRo and of OPSRo to differ from each other. 

The divergence of B2 as z^e//,r — > in eq. (|56D shows that the reduced Lifshitz equation 
(p7|) is no longer valid in this limit, i.e. as the Hopf bifurcation point is approached. At 
the Hopf bifurcation, and also at the parity-breaking bifurcations from standing structures 
to traveling structures the full complex Lifshitz equations (^) need to be used. This is not 
pursued here. 
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Additional information about possible phase-locked wave patterns is obtained by exam- 
ining the shape of the neutral surface /i^ = |z^e//P in more detail, where u^ff is given by 



eg. (^81). For example, an earlier study by Riecke, which pertains to (normal) standing rolls 
in one spatial dimension, found that the neutral curve need not be convex p]7| , p9| . In fact, 
it can have two minima, in which case the band of stable wave numbers splits into two 
disjoint regions, one centered at each of the two minima. This has been investigated, in the 
one-dimensional case, using the appropriate Ginzburg-Landau equations as well as 



phase equations [[ITipOipTI . Under the conditions that there are two minima of the neutral 
curve, it was found that wave-number gradients need not decay as would generally be ex- 
pected if there was just a single minimum. Instead, stable domain walls can arise which 
separate regions with different wave numbers q. In the present, two-dimensional study, one 
may expect that multiple minima with differing values of p can also be obtained in the 
phase-locked standing wave regime. We investigate this possibility by examining the form 
of the neutral surface, associated with the generalized Ginzburg-Landau equations (^5]), in 
a special limiting case. 

In the one-dimensional case it was shown that for multiple minima of the neutral curve to 
occur the group velocity must be small (see also equation (|5^)). Thus, to demonstrate 
the phenomena we choose f = as well as g = 0. (The minima will then persist for v 
and g sufficiently close to 0.) Figs. |^ and |] show the neutral surface /x^ = |z/e//(g,p)p for 
i>r = —0.2, d = 1 + i, s = —1 + 0.585i, and h = —4.27 + 2.5i for two different values of i>i. 
When Ui = —0.34 (fig. there is a minimum at Pc = 0, qc = 0, corresponding to normal 
PSRo, as well as a minimum (of roughly the same depth) aX p = ±Pc ~ 0.45. Thus, in 
this case one expects a competition between these three kinds of waves which could lead to 
zig-zag patterns consisting of domains with p = Pc and p = 0, respectively. In contrast to 
the zig-zags associated with the reduced Lifshitz equation (0) ||2^, one of the two zig-zag 
directions here is given by the normal direction. Changing the frequency to Ui = —0.3 shifts 
the minimum at Pc = to a non-zero value (fig. ^). Now there are four minima and one 
expects, for example, zig-zags consisting of oblique waves with two different positive values 
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of p. In the case that < 0, the extremum at g = can become a maximum and two 
new minima (in the g-direction) can arise at g = ±gi. (The reflection symmetry that forces 
symmetry-related pairs ±gi, as well as an extremum at g = 0, is due to the special choice 
V = and g = 0.) Two cases with cij < are shown in figs. ^ and ^. For di = —0.5 (fig. 
1^ the absolute minima are at g = and p = ±Pi,2! whereas for di = —0.7 (fig. ^ they are 
at p = ±pi and g = ±gi. Between these two situations the minima of the neutral surface 
almost degenerate into a circular curve around the local maximum at q = 0, p = ±0.35. 
(Recall that in the physical system g = corresponds to a non-zero wave number g.) Such a 
neutral surface, with a minimizing closed curve centered around a finite wave number, should 
lead to interesting nonlinear behavior. Of course, the complicated situations described here 
are only obtained for special parameter values {e.g. near v = 0). Generically, we expect the 
neutral surface to have only a single minimum. In particular, this is true when the group 
velocity v in (|45|) is large. 

Finally, we examine the possible stability properties of the spatially-periodic, phase- 
locked standing-wave solutions of the Lifshitz equations ( ^5D in the weakly oblique regime. 
This is accomplished by deriving the amplitude equations (^ from (^) for p ^ 0. Inserting 
the expansion 



into (|45|) , we find that the nonlinear coefficients in (^ in terms of those in (|45|) are 



with u in (^) replaced by t^eff- It then follows from our stability calculations in sections III 
and IV that PSRe are always unstable at onset in the weakly oblique regime. Moreover, 
PSRo are a stable solution at onset provided z^e//,r < and 



A = zi(t)e'(''"+P^) + Z4{t)e'^"y^ + h.o.t. 
B = 23(t)e-*(«^+™) + Z2{t)e-'^'"'-^y^ + h.o.t. 



(59) 



a = c, b = c = f = c + g, d = 2c, 



(60) 




(61) 
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This is precisely the condition that the PSRo bifurcate supercritically from the trivial solu- 
tion (c/. equation (|])). That PS Re bifurcate unstably close to the Lifshitz point is consistent 
with the observation that rectangles are an unstable solution of the reduced Lifshitz equa- 
tion ( |57D p5| , and that SRe are an unstable solution of the full Lifshitz equations ( ^5] ) in 
the absence of forcing (/x = 0) M. 



VI. COMPARISON WITH EXPERIMENTS IN ELECTROCONVECTION 

Before comparing our results with the experiments in electro-convection of nematics, it 
is useful to discuss some of the properties of the unforced experimental liquid crystal system 
that are reported in 0]. Without temporal forcing the first instability of the basic state leads 
to a steady roll pattern. Detailed investigation shows, however, strong fluctuations below 
threshold which are possibly due to thermal noise. Crucial in the present context is the fact 
that the fluctuations have the character of (oblique) traveling waves. A natural bifurcation 
diagram which is consistent with these findings is shown in fig.Q similar to convection in 
binary mixtures, there could be a subcritical Hopf bifurcation which leads to unstable waves 
which in turn merge with a branch of steady rolls. In this bifurcation scenario there is a 
discontinuous transition from the motionless state directly to finite amplitude steady rolls 
when the control parameter is increased beyond the Hopf bifurcation. 

Motivated by experiments on binary-mixture convection, where the bifurcation to trav- 
eling waves is subcritical the effect of resonant periodic forcing on a weakly subcritical 
Hopf bifurcation to one-dimensional traveling waves has been investigated theoretically . 
In this subcritical case, it was found that forcing can lead to stable PSW that arise in a 
supercritical bifurcation already below any saddle node bifurcation on the unforced (as well 
as of the forced) branch of traveling waves. In particular, bifurcation diagrams of the kind 
sketched in fig.|10|, in which the parity-breaking bifurcation of the PSW to TW is subcrit- 
ical, are possible. If the unforced traveling waves are unstable to a large-amplitude steady 
structure it is quite natural to assume that this instability persists in the case of sufficiently 
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weak forcing. Thus, at the parity-breaking bifurcation to travehng waves the system could 
jump from the PSW to the steady structure as indicated in figjlO]. 

Without temporal forcing the experimental system has two control parameters: the 
amplitude V and frequency Q of the applied ac-voltage. The bifurcation parameter z/^ 
depends linearly on the amplitude V, i.e. Ur oc {V — Vc)/Vc. The frequency Q is much 
greater than the Hopf frequency in the EHC experiments; our amplitude equations apply to 
the averaged system. The nonlinear coefficients in the amplitude equations @ depend on Q 
and thus can be varied somewhat. The resonant temporal forcing introduces two additional 
control parameters: the forcing (modulation) amplitude AV, or equivalently |yu|" oc AV'^ 
and the detuning z/j = tUe — nuh (n = 1, 2). 

In the experiments, four regimes have been identified in the (z/^, /i)-plane, as indicated in 
fig.|ll|. For small applied voltage V (above the open circles) incoherent traveling-wave patches 
are observed. They are interpreted as driven by noise, which is possibly of thermal origin 



2^ . This is substantiated by recent measurements of the spatial Fourier spectra |^2|]. Thus 
the open circles do not represent a true transition but indicate instead the location where 
the instrumental sensitivity is sufficient for the detection of very weak fluctuations. When 
the forcing amplitude is increased a transition to phase-locked standing waves is found (solid 
triangles). Depending on the driving frequency Q they are either standing oblique rolls or 
standing rectangles and can naturally be identified with the PSRo and PSRe, respectively. 
When increasing the applied voltage {i.e. u^) further, the PSRe lose stability to a structure 
(at the open diamonds) which is also phase locked to the forcing but is more complicated. 
During the phase in which the voltage rises the structure is predominantly a 'zig'. Then 
the 'zag'-component increases in strength leading to a structure which resembles rectangles 
during the peak of the voltage, whereas the 'zag'-component dominates during the decay 
phase. It is reported that this structure - in contrast to the PSRo and PSRe - shows no 
strong spatial coherence. We interpret this structure as the PSCR bifurcating off the PSRe. 
Recall that the PSCR is given by zi = = i?e**, ^2 = -^4 = re^'^, so that a typical quantity 
like the vertical velocity Vz is given by {of. (|l])) 
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Vz = 4:Rcos{uJet/n + $) cos(gx + py) + 4r cos{uJet/n + 0) cos(gx — py) + h.o.t. (62) 

Thus, provided $ — (/) is not a multiple of vr, the two standing waves differ in their phase 
with respect to the forcing; the structure alternates between 'zigs' (for values of t such that 
cos{ujf,t/n + 0) = 0), 'zags' (when cos{uJet/n + $) = 0) and rectangles (when \Rcos{uJet/n + 
$)| = \r cos{u!et/n + 0)|). Fig. |T2| shows a representative sequence of PSCR structures. 
Since the two amplitudes r and R are different, there exist two symmetry- related PSCR 
solutions, one in which the 'zigs' dominate {R > r) and one in which the 'zags' dominate 
{r > R). When the PS Re solution loses stability to PSCR, both zig- and zag-dominated 
PSCR are expected to arise and compete with each other, with domains separated by grain 
boundaries forming. It is tempting to attribute the spatial incoherence observed in the 
experiments to the existence of such domains. Unfortunately, the reported experiments do 
not provide enough data to investigate this interesting question. 

When the applied voltage is increased further a transition to a steady structure occurs 
(to the right of the solid squares) which is modulated in its amplitude. This transition 
is beyond the scope of the present analysis. One may speculate that it is related to the 
fact that the initial Hopf bifurcation seems to be weakly subcritical (c/. fig.^: without 
temporal forcing the system jumps directly to a steady branch. In the presence of forcing 
this could render the parity-breaking bifurcation of the PS Re subcritical. The system would 
then jump from the PSCR directly to the steady branch as indicated in fig.^. We have not 
studied the stabihty of the PSCR. Since the instability of the PSRe to TRe leads out of the 
standing-wave subspace, which contains also the PSCR, one might expect that they inherit 
that instability, and consequently, also exhibit a jump-transition to the steady structure. 

The present analysis suggests that the parameter p, which governs the competition be- 
tween the different standing structures, can be made small by choosing a suitable detuning 
Ui {cf. (|16|)). If this could be achieved below the unforced Hopf bifurcation (z/^ < 0), both 
transitions from PSCR to PSRo and to PSRe could be observed at small amplitudes below 
the onset of the Hopf bifurcation in the unforced case. In the experimental system the wave 
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number of the structure is, however, not fixed. Instead, the system can choose a wave num- 
ber from a continuum of values within the stabihty hmits. In particular, if v = Vr^x + Vyij is 
the group velocity in the oblique regime, then a change Agx + Apy in the wave vector shifts 
the detuning to an effective value: Veff,! = i^i+VxAq+VyAp+0{Aq'^, Ap^, AqAp). Similarly, 
near the Lifshitz point, the effective detuning at the critical wave number {QcPc = 0), with 
Qc given by (^), is strongly reduced and can remain small over a range of frequencies; in 
particular, t'e/z.i = i^i + vqc — diq^ = C(e'') [|17|. This may explain why in the experiments 
the location of the transition to PSCR, which depends on the detuning Ui or more precisely 
Veff,ii did not change substantially even though the forcing frequency was changed by almost 
50%. 

Finally, our analysis suggests that a transition to PSCR may also be found if the fre- 
quency VL is chosen such that the first transition is not to rectangles but to standing oblique 
rolls (PSRo), as long as p is positive and not too large in the regime where the PSRo arise. 
Experimentally this has not been tested in detail. 

VII. CONCLUSION 

In summary, we have studied the effect of a resonant temporal forcing on oblique traveling 
waves close to onset in an anisotropic system. We have found that such a forcing can excite 
various standing wave structures which are phase locked to the forcing: rolls, rectangles and 
cross rolls. The latter arise in secondary bifurcations from either the phase-locked roll or 
rectangle structures, and in many cases (if not all) provide a transition between them. These 
bifurcations depend strongly on external control parameters such as the forcing frequency 
(effective detuning). In particular, if the wave number of the structure could be fixed, 
one would be able to choose the control parameters such that a wide variety of bifurcation 
scenarios could be observed. In the general case, however, the attainable effective detuning is 
limited. Nevertheless, by adjusting the frequency Q of the applied ac-voltage in the unforced 
case, so that the system is close to a transition between standing rolls and rectangles, the 
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available detuning may be sufficient to observe the transition from phase-locked standing 
rolls to phase-locked standing rectangles via the general cross-rolls solution. 

A striking feature of the standing rolls and rectangles observed in the experiments is 
their strong spatial coherence. This is perhaps not surprising since they arise through bi- 
furcations in the amplitude equations that involve only real eigenvalues. Close to onset, the 
standing rolls can be modeled by a single Ginzburg-Landau equation with real coefficients. 
This equation admits a stable band of wave numbers and exhibits diffusive phase dynamics; 
moreover, it can be derived from a potential. Similarly, we expect the standing rectangles 
to be described by two coupled Ginzburg-Landau equations with real coefficients. We pro- 
pose that the spatial incoherence associated with the standing cross rolls may be due to the 
formation of grain boundaries separating domains of symmetry-related structures. A more 
detailed experimental, as well as theoretical, investigation of this regime would be interest- 
ing. Such studies could reveal whether this interpretation is correct and, if so, identify the 
dynamics associated with the domain walls. 

In the vicinity of a Lifshitz point, we model the system by two coupled complex general- 
ized Ginzburg-Landau equations which take into account the spatial degrees of freedom. We 
find, in this regime, that the phase-locked standing rectangles are always unstable. Phase- 
locked standing rolls are the preferred spatially-periodic solution in this regime. 

Finally, the results presented in this paper may also be pertinent to parametrically forced 
systems which do not exhibit a Hopf bifurcation, but which do support weakly damped 
waves. This, for instance, is the case in the Faraday experiment for low-viscosity fluids. 
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FIGURES 

FIG. 1. Phase diagram for a periodically forced supercritical Hopf bifurcation to traveling 
waves. Phase- locked standing waves bifurcate from the motionless state on the line PSW, traveling 
waves bifurcate at the line marked H. The standing waves undergo a parity-breaking bifurcation 
at PB, a saddle-node bifurcation at SN and a Hopf bifurcation to unlocked standing waves at 
SW. 

FIG. 2. Typical bifurcation diagram for a periodically forced supercritical Hopf bifurcation to 
traveling waves as obtained along the dotted line in fig.|l|. 

FIG. 3. Sketch of the bifurcation diagram of PSRo, PSRe and PSCR for A/i > and p > 
in case I. 

FIG. 4. Sketch of the bifurcation diagram of PSRo, PSRe and PSCR for A/^ < and p > 
in case I. 

FIG. 5. Neutral surface for the forced Lifshitz equation ( ^5|) with Vr = —0.2, d = 1 + i, 
s = —1 + 0.585z, h = —4.27 + 2.5i and Vi = —0.34. In the nonlinear regime zig-zag structures 
are expected to arise. The local wave number of the zigs and of the zags lies in one of the three 
minima of the neutral surface. 

FIG. 6. Neutral surface for fj = —0.3. Other parameters as in fig.^. 

FIG. 7. Neutral surface for di = —0.5. Other parameters as in fig.^. 

FIG. 8. Neutral surface for di = —0.7. Other parameters as in fig.^. 

FIG. 9. Sketch of the bifurcation diagram suggested by experiments on electro-convection in 
nematic liquid crystals. 
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FIG. 10. Effect of a periodic forcing on the subcritical Hopf bifurcation to traveling waves 
sketched in figjoj. The parity-breaking bifurcation of the PSW to TW can be subcritical, thus 
inducing a jump-transition to the steady structure. 

FIG. 11. Experimental phase diagram for temporally forced convection in nematic liquid crys- 
tals, reproduced from (See the text for a discussion of the various regions.) a) 1^ = 0.50Hz, b) 
n = 0.31Hz. 

FIG. 12. Typical time sequence for a representative phase-locked standing cross rolls solution 
(PSCR). Note that the structure alternates between a zig- and a zag-dominated state. 
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